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SYMBOLS 


A amplitude constant 

A q phasor amplitude constant 

AM amplitude modulation 

A(¥) amplitude function of approximated phase-modulated signal 

gain constant 

fundamental-frequency Fourier coefficient 
D total harmonic distortion voltage ratio 

e(t) ideal phase-modulated signal 

j imaginary 

weighting coefficients in power series 
MMSE minimum mean square error 

MTHD minimum total harmonic distortion 

n circuit order 

PM phase modulation 

Re[ ] real part of [ ] 

t time 

THD total harmonic distortion 

v(t) approximation to e(t) 

a radian frequency of modulating signal 

e 2 mean-square-error function 

<J) phase angle v(t) 

(J> o (t) phase angle of v(t) with ¥ time-varying 

m modulating variable 

design-maximum phase deviation 
peak sinusoidal phase deviation 
design-maximum peak sinusoidal phase deviation 
radian frequency of carrier signal 

iii 
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SUMMARY 


This paper examines the ultimate modulating linearity attainable with a phase 
modulation technique based on the linear addition of quadrature phase carrier sig- 
nals which have been multiplied by precisely defined nonlinear transformations of the 
modulating signal. This technique, sometimes described as complex phasor phase 
modulation, is a logical extension of the Armstrong method and can provide exception- 
ally good linearity for phase deviations as large as 5 radians. For example, it is 
shown that a total harmonic distortion ratio of -55 dB at a peak sinusoidal phase 
deviation of 5 radians is theoretically possible by using a relatively simple fifth- 
order nonlinear processor. 


INTRODUCTION 

In the early 1970's, a landing-guidance system under development required a 
phase modulator having good modulating linearity for phase deviations exceeding tt 
radians without resorting to frequency multiplication. Several methods for accom- 
plishing this objective were examined, and the one ultimately chosen is the subject 

of this paper. 

The complex phasor phase modulation technique is an extension of the Armstrong 
method (ref. 1, pp. 603 to 605) and follows logically from an expansion of the ideal 
phase modulation waveform given by equation (1) 

e (t) = Re [a 0 e j(Wt + Y) ] (1) 

into quadrature phase components. Equation (1) can be rewritten as either 

e(t) = A Q cos cot cos Y - A Q sin cot sin Y (2) 

or 

e ( t) = a (i- — + - ...) cos ut - A 0 ( V - — + — - ...) sin lot (3) 

° V 2! 4! / ° \ 3! 5! ) 

The functional configuration shown in figure 1 can be designed to perform the 
operations indicated by equations (2) and (3) . Implementation of equation (2) re- 
quires nonlinear processors to transform ¥ to sin ¥ and cos ¥. Although sine 
and cosine nonlinear function modules are available, they are two-quadrant devices, 
and their applicability is therefore limited to phase deviations not exceeding tt/2 
radians . 

In lieu of direct sine/cosine transformation, the synthesis of sine/cosine 
functions from infinite power series of ¥ is theoretically possible. This alterna- 
tive is not appealing; however, it does suggest a third possibility: the approxima- 
tion of the infinite series with a practically finite number of nonlinear terms. If 



the nonlinear processors in figure 1 can generate powers of ¥ through the fifth 
the first-zone output signal (around m) can be expressed as 


v(t) = A ^K 0 - K 2 Y 2 + K 4 ¥ 4 j cos cat - A - 


K.^ 3 + K 5 1' 5 ^ sin tot 


( 4 ) 


which may be written in the form A(¥) cos [cat + <f> 0 (¥)] with 

A OP) = A [(k 0 - K 2 v 2 + K 4 'l' 4 ) 2 + (k^ - K 3 4' 3 + K^ 5 ) 2 ] 


(5) 


and 


♦ 0 (V) = tan 


-1 


( kjY - K 3 'i ' 3 + k 5 ^\ 
K q - K 2 V 2 + K 4 Y 4 / 


( 6 ) 


This paper shows that remarkedly good linearity between $ and ¥ in 
equation (6) can be obtained if the coefficients are optimized. These "optimum" 
coefficients differ significantly from those in equation (3) . This paper addresses 
the following questions: 

(1) What optimization criteria should be used? 

(2) What are the "optimum" coefficients in equation (6 ) , based on the criteria 
established in (1) above? 

(3) What total harmonic distortion performance can be expected when the optimum 
coefficients are used? 

Since the initial work on this project was done over a decade ago (ref. 2) , 
significant improvements in the components needed to implement the radio-frequency 
and nonlinear processors shown in figure 1 have occurred (ref. 3 and ref. 4, p. 13) , 
and further improvements are likely. Therefore, this study is limited to defining 
the ultimate performance capabilities of the technique without regard to hardware 
imperfections. Circuit orders of 2 through 5 will be considered, where order denotes 
the highest power of Y having a nonzero coefficient in equation (6) . 

The fact that the envelope of the synthesized approximation to e(t) is not 
constant but is a function of Y as given by equation (5) is largely ignored, and it 
is assumed throughout most of the paper that hard limiting can satisfactorily remove 
all amplitude fluctuations. A final section describes the incidental envelope varia- 
tions that inevitably accompany this phase modulation technique for the purpose of 
facilitating the design of a suitable limiter. 


ANALYSIS 

Phase Function Optimization 

The approach employed here was to postulate reasonable and practically 
applicable optimization criteria, determine the optimum Kj_ set based on each of 
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these criteria, and determine the linearity between <f> Q and T for each set on 

a common basis. The optimization criteria utilized here fall into the general cate- 
gories of minimum mean square error (MMSE) and minimum total harmonic distortion 
(MTHD) for T, a sinusoidal signal. The final choice between the various approxima- 
tions was based on total harmonic distortion (THD) performance. That is, for T (t) 

= AT sin at, the Fourier coefficients of cf) Q (t) were computed as a function of 
AY for each set. It should be emphasized that THD minimization to determine 

the optimum Kj_ set and the THD comparison between the various approximations (Kj_ 
sets) are quite different operations. In the first case, AT is fixed at AT D and 
the THD is minimized with respect to the values, whereas in the second case, the 

K-j_ values are fixed and AT is the independent variable. 

The mean-square-error (MSE) function is defined as 



where C Q T is the desired linear relationship between <j> 0 and T. In general, the 
mean-square-error function is minimized with respect to the values; however, 

other constraints or conditions may be imposed. An initial optimization approach, 
designated (MMSE) set Kq and in equation (6) equal to unity, and the coef- 

ficient C Q in equation (7) was made equal to the first derivative of <}> 0 (T) in 
equation (6) evaluated at 0, or unity. Thus, for an nth-order circuit, the procedure 
determines the n - 2 remaining coefficients. Computer programs were developed which 
minimized the integral expression in equation (7) by a numerical procedure related 
to the Newton-Raphson method (ref. 5, p. 221) generalized to an arbitrary number of 
independent variables. The integrals were evaluated numerically by means of 
Simpson* s rule. In figures 2 through 5 the optimum (MMSE) g coefficients for circuit 
orders 2 through 5 are plotted as a function of T^, the design-maximum phase 
deviation. 

The normalized root mean square error, or RMSE , provides a physically meaningful 
measure of the deviation of 4> Q (T) from the desired linear relationship C 0 T. 
Normalized RMSE is defined as 

( 8 ) 

2 ^ 

and is plotted versus T D in figure 6 for the purpose of indicating the circuit 
order required to achieve a specified value of normalized RMSE for values of T D 
up to 5 radians. 

Total harmonic distortion ratio, a widely accepted measure of system 
nonlinearity, is defined as 


Normalized RMSE = 


A 

D = 


Total harmonic distortion "power 11 
Fundamental "power" 


( 9 ) 
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when the excitation is sinusoidal, that is, when T(t) = AT sin at. The THD ratio 
may be expressed in terms of the mean square value of (p Q and C-^, the fundamental- 
frequency Fourier coefficient, as 




(10) 


If 'F (t) is sinusoidal with period 2n/a, <f> 0 (t) also has a period of 2ir/a and 
possesses both odd and half-wave symmetry. The total "power" in <(> ( t) is 


-ir/a 

> 2 = — f 

'o 2 IT J 


<t> 0 (t) dt 


(ID 


-n/a 


and the fundamental- frequency "power" is 




n2 


— -ir/a 

J ^o (t) sin at dt 


-■n/a 


(12) 


An integral expression for D 2 (A’l') is obtained upon combining equations ( 6 ), 
(10), (11), and (12), and letting 'l'(t) = A'F sin at. This expression was evaluated 
numerically by using the (MMSE) optimum sets corresponding to selected values 

of Tp. The results are plotted in figures 7 through 10 for circuit orders of 2 
through 5. The THD performance of the Armstrong (linear) approach is shown as a 
dashed line for comparison purposes. The THD was plotted up to values of AT 
numerically equal to the Tp values producing approximately 1 percent normalized 
MSE in figure 6 . The following observations can be made about these curves: 

(1) The THD performance is much improved relative to the Armstrong method. 

(2) In all cases, THD increases very rapidly when AT exceeds Tp. 

(3) As AT becomes less than Tp, THD fluctuates generally downward; however, 
THD is always less than the value corresponding to AT = Tp. 

(4) When AT « T D , THD decreases monotonically with AT. 

(5) Values of THD as low as -55 dB can be obtained with a fifth-order circuit 
for a peak sinusoidal phase deviation of 5 radians, as shown in figure 10. 

In addition to (MMSE) ^ minimization, several other optimization criteria were 
postulated and evaluated but were rejected in favor of (MMSE)i- Additional MMSE 
constraints were imposed in an attempt to further reduce the normalized MSE. In 
(MMSE ) 2 minimization, in equation ( 6 ) was made an additional independent varia- 

ble and equation (7) was minimized, with 0 o = 1. For (MMSE ) 3 minimization, C Q in 
equation (7) was made the additional independent variable and equation (7) was mini- 
mized, with = 1. In both cases, the reduction in normalized RMSE was trivial; 
however, the THD behavior was considerably different from that of (MMSE)i, as 
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discussed shortly. In THD minimization, the integral expression for D 2 (AY D , Kj_) , 
obtained by combining equations ( 10 ) , ( 11 ) , and ( 12 ) , was minimized with respect 
to the Kj_ values, with AY D being a parameter. Thus, D 2 was minimized for 
Y(t) = A¥ d sin at. Upon obtaining the optimum Kj_ set for each value of A¥ D , 
D 2 (AY) was computed as in the (MMSE) ^ case. 

The salient results of comparing D 2 (AY) for all four optimization criteria 
can be summarized as follows (refer to fig. 11 ) : 

(1) For AY « AY D (or Y D ) , (MMSE) ]_ is slightly superior to the other three 
cases. 

( 2 ) For intermediate values of AY, say AY = 1/2AY D , (MMSE) 2 , (MMSE) 3 , and 
MTHD all exhibit an undesirable characteristic - THD increases considerably 
relative to the THD value corresponding to AY = Y D (or AY D ) . 

(3) For values of AY approaching AY D , MTHD always gives the lowest THD, 

as it should; however, the improvement is small - typically less than 5 dB. 

Figure 11 illustrates these effects and is representative of the behavior 
observed with other circuit orders and values of AY and Y D . Thus, the small re- 
duction in THD at values of AY approaching Y D ( (MMSE ) 2 and (MMSE) 3 ) and AY D 
(MTHD) was more than offset by the undesirable behavior at intermediate values of 
AY and the larger THD at small values of AY. Therefore, (MMSE)]_, the simplest 
both conceptually and computationally, was judged to be the best choice for 
calculating the weighting coefficients in equation (4) . 


Concomitant Envelope Variation 

The envelope function of equation (5) was totally disregarded in linearizing 
the phase function of equation ( 6 ) . Thus, it should not be surprising to find that 
appreciable envelope variation accompanies this phase modulation method. Well- 
designed hard limiting, as indicated in figure 1 , can greatly reduce residual enve- 
lope variations; however, AM to PM conversion in the limiter (ref. 6 , p. 51) should 
be strictly avoided. 

The degree of envelope variation to be encountered with this PM method is 
therefore of interest. This information is contained in figure 12, which shows 
A(Y)/A(0) (see eq. (5)) as a function of the independent variable Y for various 
circuit orders and values of Yp. It is seen that A(Y) is a monotonically increas- 
ing function of Y over the range of variables considered, and that higher circuit 
orders always result in less envelope variation. 


EXPERIMENTAL RESULTS 

As noted in the Introduction, this phase modulation technique was applied to a 
practical situation. A fifth-order circuit for Y D = 5 radians was built, evaluat- 
ed, and successfully used in the intended application. Unfortunately, the apparatus 
available for evaluating the modulating linearity of the experimental phase modula- 
tor did not possess adequate linearity to make a definitive measurement. The theory 
predicted that the experimental model should exhibit THD not exceeding -55 dB (see 
fig. 10 ) ; however, measurements made with the test apparatus available resulted in 
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THD values of -40 dB . This discrepancy is at least partially, and perhaps largely, 
the result of distortion generated within the phase demodulator used to evaluate the 
experimental phase modulator. Thus, the extent to which a practical ’circuit matches 
the ideal remains undefined and may well be better than the measured value of 
1 percent, or -40 dB. Even if the -40-dB value is correct, it is the opinion of the 
author that the hardware improvements mentioned in the Introduction would result in 
better THD performance. 


CONCLUDING REMARKS 

This paper has shown that the complex phasor phase modulation technique is 
capable of exceptionally good modulation linearity for phase deviations as large as 
5 radians. A simple procedure was developed to determine the optimum weighting 
or gain coefficients in the nonlinear processors required to implement the technique. 
These optimum coefficients are shown graphically for circuit orders up to 5 and phase 
deviations up to 5 radians. The total harmonic distortion performance for a complex 
phasor modulator using the optimum coefficients is defined. Finally, the unavoidable 
amplitude modulation accompanying the technique is quantified. 

Although this paper is basically theoretical in nature, the ideas contained 
herein were applied, and practical circuits were constructed. For those familiar 
with analog and radio-frequency signal-processing components, it will be apparent 
that the critical items in this concept are (1) the analog multipliers needed to 
implement the nonlinear processors and (2) the balanced modulators (also analog 
multipliers) used to perform the carrier modulation operations. Since the other 
hardware items needed to implement the concept can be considered to be essentially 
perfect for this application, the ultimate performance will be dependent on the 
analog modulators. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
June 24, 1985 
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Figure 1.- Functional diagram of complex phasor phase modulator. 
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Figure 2.- (MMSE) ^ K 2 values. Circuit orders 2 through 5. 
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Figure 3.- (MMSE) values. Circuit orders 3 through 5. 
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Figure 5.- (MMSE) ^ values* Circuit order 5. 
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Figure 7.- Total harmonic distortion versus peak phase deviation 
Second-order circuit; (MMSE) n coefficients. 
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Figure 9.- Total harmonic distortion versus peak phase deviation 
Fourth-order circuit; (MMSE) , coefficients. 
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Figure 10,- Total harmonic distortion versus peak phase 
Fifth-order circuit; (MMSE) ^ coefficients. 
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Figure 11.- Total harmonic distortion for various optimization 
criteria versus peak phase deviation. Third-order circuit. 
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